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Cosmological variation of the fine structure constant a due to the 
evolution of a spatially homogeneous ultra- light scalar field (m ~ Hq) 
during the matter and A dominated eras is analyzed. Agreement of 
Aa/a with the value suggested by recent observations of quasar ab- 
$_i ■ sorption lines is obtained by adjusting a single parameter, the coupling 

c/2 . of the scalar field to matter. 

Asymptotically a(t) in this model goes to a constant value a ~ ao 
in the early radiation and the late A dominated eras. The coupling of 
the scalar field to (nonrelativistic) matter drives a slightly away from 
a in the epochs when the density of matter is important. 

Simultaneous agreement with the more restrictive bounds on the 
variation |Aa/a| from the Oklo natural fission reactor and from me- 
teorite samples can be achieved if the mass of the scalar field is on the 
order of 0.5-0.6 Ha, where H\ = Q^ 2 Hq. 

Depending on the scalar field mass, a may be slightly smaller or 
larger than ao at the times of big bang nucleosynthesis, the emission of 
the cosmic microwave background, the formation of early solar system 
meteorites, and the Oklo reactor. The effects on the evolution of a 
due to nonzero mass for the scalar field are emphasized. 

An order of magnitude improvement in the laboratory technique 
could lead to a detection of (a/a)o. 



1 Introduction 



Recent observations by Webb et al. [TJ|2] of absorption lines in quasar spectra 
provide evidence for a variation of the fine structure constant 

Aa a(t) — an , = 

= — ^ = -0.57 ± 0.10 x 10"° 1 

a a 

averaged over the redshift range 0.2 < z < 3.7 ("a was smaller in the past"), 
where a is the present-day value of the fine structure constant. This type of 
variation of a, as well as variation of other dimensionless coupling constants, 
is predicted by theories which unify gravity and the standard model forces. 
For example, string and supergravity theories predict the existence of mass- 
less or ultra-light scalar fields (dilaton or moduli fields) which through their 
dynamical evolution can cause temporal variation of coupling constants. 

This investigation will consider cosmological variation of the fine structure 
constant due to the evolution of a spatially homogeneous ultra-light scalar 
field (m ~ H , where H is the present value of the Hubble parameter) 
during the matter and A dominated eras. We will assume a flat Friedmann- 
Robertson- Walker universe, with p c0 = p m o+p r o+PA ~ Pmo+PA today, where 
p c o is the present value of the critical density for a flat universe, and p m n, 
Pro <C Pmo , and p\ are the present energy densities in (nonrelativistic) mat- 
ter, radiation, and the cosmological constant respectively. Ratios of present 
energy densities to the present critical density are denoted by Q m = p m o/Pco, 

= Pro/Pco, and fi A = pa/pco- 

The scalar field (ft may provide the cosmological constant energy density at 
the minimum of its potential V((ft). In the model presented here, the energy 
density p^ — V(<f>) of the scalar field is always very small compared with the 
critical energy density in the radiation, matter, and A (dominated) eras, so 
that the standard Friedmann- Robert son- Walker evolution of the universe is 
not affected by displacements of (ft from (ft. 

Agreement of Aa/a with the quasar data can be obtained by adjusting a 
single parameter, the coupling of the scalar field to (nonrelativistic) matter. 
Asymptotically a in this model goes to a constant value a ~ ao in the early 
radiation and the late A eras, insuring agreement with bounds from cosmic 
microwave background (CMB) temperature fluctuations (| Act/a | < 0.05 at 
z = 1090) and big-bang nucleosynthesis (BBN) (\Aa/a\ < 0.02 at z ~ 10 9 - 
10 10 ) (see Ref. |3j for a comprehensive review). 
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Simultaneous agreement with the more restrictive bounds on the total 
change |Aa/a| < 10~ 7 from z m 0.14 to the present from the Oklo natural 
fission reactor 1.8 Gyrs ago |3J Ej (by analyzing isotopic ratios of Sm) and 
|Aa/a| < 3 x 10~ 7 from z ~ 0.44 to the present from samples of meteorites 
formed in the early solar system 4.6 Gyrs ago jH] (by analyzing the ratio of 
187 Re to 187 0s) can be achieved if the mass of the scalar field is on the order 
of 0.5-0.6 n]! 2 H . 

The laboratory bounds on the present variation |d/a|o < 3.7xl0~ 14 /yr [7] 
and (d/a)o = (—0.4 ± 16) x 10~ 16 /yr t SJ are satisfied for the entire range of 

1 /2 

scalar field masses < m < 12 fi A Hq considered here. The variation (d/ajo 
predicted in the model may be detectable if the sensitivity of the laboratory 
experiments can be increased by an order of magnitude. 



2 The Model 

The scalar field model is based on a generalization of Bekenstein's model 
for variable a, but with an ultra-light scalar field mass. The scalar field obeys 
the evolution equation 

in the standard Friedmann- Robertson- Walker cosmology. Here H is the Hub- 
ble parameter, ( m is the coupling of to matter, \( m \ <C 1, p m is the den- 
sity of matter, M* $ Mp is the mass scale associated with the scalar field, 
and the (reduced) Planck mass Mp = 2.4 x 10 18 GeV. General considera- 
tions show that the coupling of <fi to radiation (including relativistic matter) 
should vanish — since (f) couples to the trace of the energy-momentum tensor 
for matter and radiation — and that ( m is very nearly constant during the 
matter and A eras. 

In generalizations of Bekenstein's model, variation of a derives 

from the coupling of to the electromagnetic field tensor F^, through a 
term in the action of the form 

S F = J d A xy^ (-^Bp^/M^F^ (3) 

where Bp is a function (introduced by Damour and Polyakov [TH]) that would 
be specified by the string or supergravity theory and constitutes the effective 
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vacuum dielectric permittivity. In Bekenstein's model, Bp can be written as 
Bp = exp{— 2(0 — 0)/M*}. Changes in <p induce changes in a: 



a(t) 



a 



B F ((f>(t)/M*) 



with B F {<t)/M*) = 1. 

Our attention will be restricted to small departures of 
will occur in the radiation, matter, and A eras. Defining 



(4) 



from d> which 



<P 



the equation for the evolution of the scalar field becomes 



if + 3Hip + m 2 (p = —( Ti 



Pm 



j- Pm 



(5) 



(6) 



to first order in ip, where m 2 = V" ((/>), ( = Mp( m /M 2 , \(\ < 1, a is the scale 
factor, and ao is its present value. For small (p, Eq. (J2J) becomes 



a(t) w a (1 - (fv) 



and 



Aa 

~ Cf(<Po 

a 



(7) 



(8) 



where (p = B' F ((f)/M*). In Bekenstein's theory, ( F = —2. 

The experimental constraints from the validation of the weak equivalence 
principle on the couplings ( and (p may be evaded by assuming that couples 
predominantly to dark matter ^J^j. (For a different view, see the extended 
discussion of varying a and the equivalence principle tests in Ref. [To].) 

Given a complete particle theory, (p will be specified and it will be possi- 
ble to calculate the coupling £ of to matter. However, the sign and magni- 
tude of ( vary depending on the way in which Bekenstein's theory is gener- 
alized ^2] — and can depend on the unknown properties of dark matter — so 
here ( will simply be determined to fit the quasar data. 

One way in which an ultra-light scalar field mass might arise is that near 
de Sitter space extrema in four-dimensional extended gauged supergravity 
theories (with noncompact internal spaces), there exist scalar fields with 
quantized mass squared 



m 



nHl Hi 



Pa 
3Ml 



(9) 



4 



where — 6 < n < 12 is an integer and if a is the asymptotic de Sitter space 
value of H with cosmological constant p^. In certain cases, these theories are 
directly related to M/string theory. An additional advantage of these theories 
is that the classical values m 2 = nH\ and p\ are protected against quantum 
corrections. (Cosmological consequences of such ultra-light scalars in terms 
of the cosmological constant and the fate of the universe are discussed in 
Refs. HH1 1221 123) 

Note that the relation m 2 = nH\ was derived for supergravity with scalar 
fields; in the presence of other matter fields, the relation may be modified. 

We will take n > 0, corresponding to a de Sitter space minimum, and 
will contrast the evolution of a with n > with the massless case n = 0. 
For n > 0, ip — > and consequently a — > a as t — > oo. It is always possible 
to satisfy the quasar constraints on Aa/a for integer < n < 12, except 
for n — 1. However, the limits on variation of a from the analyses of Oklo 
and meteorite data are not simultaneously satisfied in this model unless n = 
0.24-0.34 (m w 0.5-0.6 H A ). 

A supergravity inspired potential [213 E] for the scalar field is 

Vffl = Pa cosh ( ^J^ ) = PA cosh(y,) (10) 

with M* = Mpj\pl. This potential produces the present-day cosmological 
constant p\ when m <p and an ultra-light scalar field mass 



- 2 (0 = 0) = !^ = 6tfl. cm 



The potential (fTU|) provides a specific realization of the generic case (0), with 
V((p) = pa- Even in this case V{<p) may have a more complicated form in 
general and only approach pACOsh(<p) asymptotically, for example, after a 
symmetry breaking phase transition. (For an analysis of varying a in models 
with a "quintessence" potential V((p), see Refs. (2HI2S1-) 

A major difference between the present model and that of Refs. |H] - P2! is 
that here the scalar field is assumed to be near the minimum of its potential, 
and thus a(t) — > a for t Hq . The initial conditions advocated below also 
differ from those of Refs. |H]-[I2], and insure that a always remains close to 
a. 

While a mass term is allowed in the generalized model of Olive and 
Pospelov ^2], h is neglected for the explicit solution given there in Eq. (3.4). 
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In a later section, the authors consider a mass term in the context of the 
Damour-Polyakov model ^3] for varying constants, where ( m = = (p. In 
this model, 

^ - \& - (12) 
where £p = B F ((f)/M*), and the scalar field mass is given by 




The authors find that the variation of a can be made marginally consistent 
with the quasar and Oklo data if the £a term in m 2 is positive and dominant 
over the £ m term. 

By contrast, the detailed effects on the evolution of a due to nonzero 
mass for the scalar field with ( F ^ and ( m ^ will be emphasized below. 

3 Evolution of the Scalar Field 

To determine the initial conditions for the evolution of the scalar field, we 
will match approximate solutions to the evolution equation © from the 
radiation era and the matter era at the time t m _ r (z ~ 3200) of matter- 
radiation equality. Note that in the early radiation era, the right-hand side 
of the evolution equation (JBJ) goes to zero, since ( m = for radiation. 

In the early radiation era, <f> may be displaced from cf) and "frozen" due to 
the large frictional term 3H<fi in the evolution equation. However, for nonzero 
m ~ H , the magnitude of the initial value \(pi\ in the early radiation era must 
still be 1 to satisfy the BBN, CMB, and quasar bounds on the variation 
of a. To see this, note that for m ~ H Q , tp is frozen near (ft until H becomes 
of order H and then decays with a characteristic timescale on the order of 
Hq . Thus the change in from the early radiation era to the present and 
the concomitant change Aa/a are on the order of tpi (see Fig. ^). To satisfy 
the quasar bounds on Aa/a would require fine tuning of initial conditions 
(^^-10" 5 ). 

For m — 0, the initial condition for (p is irrelevant in the linearized the- 
ory © since Act only depends on changes in <p; the value of ipi becomes 
important only if the scalar potential cannot be neglected. 

It is plausible though that if = in the early radiation era, since V(<j)) may 
have a deep minimum at during inflation in the very early universe, which 
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later, after one or more phase transitions, becomes shallow with V"{$) ~ Hq. 
For example, in the primordial inflationary stage, V"(4>) may be on the order 
of H], where Hi 3> Hq is the Hubble parameter during primordial inflation. 
While the scale factor a inflates by 60 or more e-foldings, will in this 
scenario rapidly approach <p. We therefore take cp — = (p in the early 
radiation era. 

Equation (jOJ) may be put into dimensionless form by setting r = H\t = 
^ + 3^ + ^ = -3C^(^) 3 (14) 

JJA "A V 0, J 

where henceforth a dot over (p denotes differentiation with respect to r. 

The Hubble parameter is related to the scale factor and the energy den- 
sities in matter, radiation, and the cosmological constant through the Fried- 
mann equation 



H\ \adtr J Vt^ \ a ) ' Qa V a 
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which can be used to solve for a and H in the matter-A, early matter, and 
radiation eras. 

In the matter-A era, the scale factor and Hubble parameter have the 
explicit forms 



t ^ sinh 2 / 3 -r (16) 

a i + z \n A J 1 1 1 

H = # A coth f-r) (17) 



and the evolution equation becomes 



'3 \ —3( 

(p + 3 coth ( -r K? + nip = 9 /o \ • ( 18 ) 

/ sinh 2 (Jr) 



In the early matter era, the scale factor 

/^A 1 / 3 /3r\ 2 / 3 



(19) 



and the Hubble parameter H/H\ « 2/(3r). The mass term nyj in Eq. 
can be neglected in the early matter (and radiation) eras. The evolution 
equation for the scalar field in the early matter era becomes 

2 AC 

V + ~V = -^- 2 (20) 

which has the solution 

( Pm = ~ (in — + - + <*) (21) 

"J \ 'Tm—r T / 

where r m _ r = H/^t m _ r and ci and c 2 are constants of integration. 
In the radiation era, the scale factor 

AQ \ 1/4 , 



a V 



s = T _^zL (22) 



and the Hubble parameter H/H^ « l/(2s). The time shift r m _ r /4 is de- 
termined by matching the Hubble parameters from the radiation and early 
matter eras at r m _ r . (The formal mathematical singularity at a — > now 
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occurs at t = t m _ r /A due to the choice of the zero of time in Eq. (JTBJ) .) The 
evolution equation for the scalar field becomes 

Cp + — <p = ~ 3C . (23) 
V 2s 4s 1 / 2 s 3/2 v ; 

The solution in the radiation era is 

<Pr = ~¥- I + (24) 

^m—r 

where C3 and C4 are constants of integration. 

The initial conditions for the scalar field in the early radiation era are 
(friji) = = <£v( r t)> where the initial time 73 satisfies s» <C s m _ r . These 
initial conditions fix the constants in the solution (|24|) . yielding 

* 771 — 7* TYl — 7" ^ m — 1 



Next take the limit Sj/s m _ r — ► to obtain 

3C s 1/2 

Vr « -y ITT • ( 26 ) 

Now match </v = y9 m and ip r = ip m at r m _ r to determine the constants in 
the solution (JUJ): 

*- = -ffob + 2 Zr + I) < 27 > 

4C/1 r m _ r \ 

^ = -~^T • ( 28 ) 



3 Vr 4r 2 

To simulate the evolution of the scalar field in the matter-A era, we use 
Eq. (fTS|) with initial conditions provided by Eqs. (J2~7j) and (|25|) evaluated at 

Tin. — t • 
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4 Comparison with Quasar, Meteorite, and 
Oklo Data 



In numerical values for expressions, we take Vt m = 0.27, Q\ = 0.73, and H 
= 71 (km/sec) /Mparsec = 1.5 x 10~ 33 eV, from Table 10 of the first-year 
WMAP observations With these parameters, the age of the universe is 
to = 13.7 Gyrs, and the absorption clouds at z = 0.2-3.7 date to 2.4-11.9 
Gyrs ago. 

Figs. I2TH31 present simulations of the evolution of the scalar field and 
Aa/a for n = 6, 12, 2, 1, 0, and 0.3. 

In the figures for Aa(z)/a, the dark (1 a error bounds for 1 < z < 2.5) 
and light (a rough guide to the error bars for 0.6 < z < 3) boxes indicate the 
quasar bounds from the bottom panel of Fig. 2 of Ref. [2], while the short 
vertical lines at z = 0.14 and z = 0.44/2 indicate the Oklo and approximate 
meteorite bounds respectively. While the Oklo event is sharply located in 
time, the ratio of 187 Re to 187 0s observed today in meteorites involves the 
total change in a since the time of formation of the meteorites, which is 
calculated by averaging a — a®. 



The scalar field (p solving the initial value problem defined in Eqs. (|18|) . 
(J2ZJ), and (|2~%j) will be proportional to (, and thus Aa/a will be proportional 
to ((f- For simplicity we will set (f = — 2, but a general value for (p 
can be reinserted. For < n < 12 except n = 1, (is fixed by setting 
Aa/a = —0.57 x 10 -5 at z — 1.75. For n = 1, better results were obtained 
by setting Aa/a = —0.57 x 10~ 5 at z = 1 (see Fig. EJ). The massless case 
shown in Fig. ITU1 agrees with Eq. (3.4) of Ref. [12] with £a = 0. For n = 0.3, 
Figs. El and El show that the quasar, meteorite, and Oklo bounds can be 
satisfied simultaneously. 

The number of visible oscillations in the scalar field p (and thus also in 
Aa/a) corresponds to how massive the scalar field is, with at one extreme 
no oscillations for the massless case fFig. I10J). and at the other extreme two 
visible oscillations for the n = 12 case (Fig. 

In this model, for n > 1 (n = 1), a was actually larger in the past at some 
point before the period of the z = 0.2-3.7 (z = 0.2-1.9) absorption clouds, 
and will be larger again in the future. For n = 0, a was smaller in the past, 




(29) 
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but will be very slightly larger in the future. And for n = 0.3, a was smaller 
in the past until just after t Q kio, was slightly larger from that point up to t , 
and then again will be smaller in the future (Aa/a — > — 10~ 4 ). 

Values of £, BBN, CMB, meteorite, and Oklo Aa/a, and (d/a)o for 
various scalar field masses are presented in Table 1. 

The BBN, CMB, and quasar bounds on Aa/a and the laboratory bound 
on |d/a|o are satisfied for integer < n < 12, except that the n = 1 case 
cannot be made to satisfy the quasar bounds in this model. The variation 
\Aa/a\ satisfies in addition the Oklo bound for 0.26 < n < 0.34 and the 
meteorite bound for 0.24 < n < 0.43 (see Fig. HI)). There is a small range of 
scalar field masses for which \Aa/a\okio < 10~ 8 (and can be made to go to 
zero by extreme fine tuning). For this range of n, (Aa/a) meieor ~ — 1 x 10 -7 . 

Depending on the scalar field mass, the predicted BBN, CMB, meteorite, 
and Oklo values of a may be slightly smaller or larger than ao- Note that 
the sign of ( for n = is opposite to the sign for n > 1. An order of mag- 
nitude improvement in the experimental technique could lead to a detection 
of (a/a) . 

For the massless case, the variation in a can be made marginally con- 
sistent with the quasar, meteorite, and Oklo bounds by setting Aa/a = 
-0.18 x 10~ 5 at z = 3 (Fig. HU). 

The behavior of a (z) can pin down the values for ( and m. Conversely, 
even knowing only the sign of ( or Aa/a can rule out certain values of the 
scalar field mass. For example, if (A«/a) BBJ v > or (( F < 0, then n ~ 1. 

5 Conclusion 

Asymptotically a(t) in this model goes to a constant value a ~ ao in the 
early radiation and the late A dominated eras. The coupling of the scalar 
field to (nonrelativistic) matter drives a slightly away from a in the epochs 
when the density of matter is important. 

Even for = 0,0:^57 as t— > oo as long asm^O. In the massless 
case, as t — ► oo, a goes to a constant value which differs from a but still 
approximately equals 1/137 for JXv > 0, while if fl\ = (and V(4>) = 0), 
Aa ~ (Xi?hi(T/To), as in Ref. jTTj. Thus the variation |Aa/a| of the fine 
structure constant becomes of order 1 only if both p\ — ■> and m = 0, and 
only for t ^> t . 

The simulations above indicate that it is possible to extract properties of 
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the scalar field from quasar absorption line spectra, including the coupling 
of (j) to matter and its mass. The variation of a has different behaviors in 
the redshift range < z < 5 depending on the mass of the scalar field. 
Thus additional quasar absorption line data, and better Oklo and meteorite 
bounds, will help elucidate the properties of the scalar field. The case m = 
H\ is ruled out in this model. A laboratory detection of (a/a) may be 
possible in the near future. 

To satisfy the quasar, meteorite, and Oklo bounds on Aa/a, the mass of 
the scalar field has to be on the order of 0.5-0.6 H^. It is difficult to satisfy 
both the Oklo/meteorite and quasar bounds in theories where the variation 
of a derives from the evolution of a scalar field; the scalar field in the model 
studied here must be near an extremum near t Q kio and t . 

The key insight of this model, as well as other models of variable a, is 
that variation of a provides a window into the parameters of the underlying 
theory that unifies gravity and the standard model of particle physics. 
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n 
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\'— iu -/ u -)rneteor 
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(a/ a)o 


n 

u 


9 n v i n~ 6 

— Z.U A 1U 


— U. i A 1U 


— O.O A 1U 


— U.O A 1U 


—A R y 1 n~ 7 


9 A y 1 n~ 16 


0.3 


-3.0 x 10~ 6 


-9.6 x 10~ 5 


-7.6 x 10~ 5 


-1.4 x 10~ 7 


-3.2 x 10~ 9 


-6.8 x 10~ 17 


0.31 


-3.0 x 10~ 6 


-9.8 x 10~ 5 


-7.7 x 10~ 5 


-1.1 x 10~ 7 


2.0 x 10~ 8 


-8.3 x 10~ 17 


1 


6.1 x 10~ 6 


1.8 x 10~ 4 


1.4 x lO" 4 


-3.7 x 10- 6 


-3.2 x lO" 6 


2.0 x 10~ 15 


2 


1.9 x 10~ 6 


5.0 x 10~ 5 


3.8 x 10~ 5 


-2.7 x 10~ 6 


-2.2 x 10~ 6 


1.3 x 10~ 15 


3 


1.0 x 10" 6 


2.4 x 10~ 5 


1.7 x 10" 5 


-2.2 x 10- 6 


-1.8 x 10- 6 


1.1 x 10~ 15 


4 


7.0 x 10~ 7 


1.5 x 10" 5 


1.0 x 10" 5 


-2.0 x 10~ 6 


-1.6 x 10~ 6 


9.3 x 10~ 16 


6 


4.6 x 10~ 7 


7.6 x 10- 6 


4.4 x 10~ 6 


-1.8 x lO" 6 


-1.4 x 10~ 6 


8.0 x 10~ 16 


12 


2.6 x 10~ 7 


1.5 x lO" 6 


-2.8 x 10~ 7 


-1.4 x 10~ 6 


-1.1 x 10~ 6 


5.7 x 10~ 16 



Table 1: Values of (, BBN, CMB, meteorite, and Oklo Aa/a, and (a/a)o in 
yr _1 vs. scalar field mass squared n. 
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Figure 2: Scalar field ip and Aa/a vs. £/£q for n = 6. 
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Figure 3: Aa/a vs. z for n = 6. 
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Figure 9: Aa/a vs. z for n—1. 
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Figure 10: Scalar field ip and Aa/a vs. t/to for n = 0. 
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Figure 11: Aa/a vs. z for n — 0. 
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Figure 12: Scalar field ip and Aa/a vs. t/to for n = 0.3. 
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Figure 13: Aa/a vs. z for n = 0.3. 
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Figure 14: Aa/a vs. n satisfying the Oklo and meteorite bounds. 
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Figure 15: Aa/a vs. z for n = by setting Aa/a = —0.18 x 10~ 5 at z — 3. 
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